Limits of the equivalence of time and ensemble averages in shear flows. 
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In equilibrium systems, time and ensemble averages of physical quantities are equivalent due 
to ergodic exploration of phase space. In driven systems, it is unknown if a similar equivalence 
of time and ensemble averages exists. We explore effective limits of such convergence in a sheared 
bubble raft using averages of the bubble velocities. In independent experiments, averaging over time 
leads to well converged velocity profiles. However, the time-averages from independent experiments 
result in distinct velocity averages. Ensemble averages are approximated by randomly selecting 
bubble velocities from independent experiments. Increasingly better approximations of ensemble 
averages converge toward a unique velocity profile. Therefore, the experiments establish that in 
practical realizations of non-equilibrium systems, temporal averaging and ensemble averaging can 
yield convergent (stationary) but distinct distributions. 
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A central tenant of statistical mechanics is the ergodic 
hypothesis: during time evolution, thermodynamic sys- 
tems pass through almost all possible microstates lead- 
ingto the equivalence of time and ensemble averages 
(Hill IE 0> Hj] • The exploration of phase space is achieved 
through thermal fluctuations, and this allows for a de- 
scription of equilibrium dynamics in terms of well defined 
averaged observables. The observables are extracted ei- 
ther by averaging over states explored over long periods 
of time or from ensembles of states distributed over all of 
phase space, whichever method is most convenient for the 
observable being studied. A natural questions is whether 
or not athermal fluctuations in a driven system lead to 
the equivalence of time and ensemble averages through 
a similar exploration of phase space. Ideally, testing the 
equivalence of time and ensemble averages in driven sys- 
tems requires measuring the instantaneous value of an 
observable for an infinite period of time (temporal aver- 
ages) and constructing an infinite number of ensembles 
(ensemble averages). This is untenable for real experi- 
ments; however, if time scales of observation are faster 
than intrinsic timescales, it is possible to gather a rea- 
sonable approximation of the distribution of values for 
an observable. 

The averaging process is governed by the fluctuations 
that move the system through phase space. In driven, 
complex fluids, athermal fluctuations arise through par- 
ticle rearrangements that are induced by flow. The po- 
tential for these athermal fluctuations to play a role 
that is similar to thermal fluctuations in normal fluids 
is the basis of many statistical treatments of highly non- 
equilibrium complex fluids A specific example 
of this are the various proposals for using effective tem- 
peratures @, 0, H| in driven systems. An open question 
in this field is the quantitative disagreement [lfj between 
different definitions of effective temperatures [9(. These 
disagreements may be connected to the ability of ather- 
mal fluctuations to provide sufficient exploration of phase 



space for certain variables, such as velocity. In this letter, 
we specifically focus on the ability of athermal fluctua- 
tions in the system to explore phase space in a manner 
that results in an equivalence between time and ensem- 
ble averages. Testing such an equivalence is an important 
step in understanding how well the various elements of 
statistical mechanics (such as temperature) can be trans- 
lated to driven systems. 

In considering any comparison of time and ensemble 
averages, it is important to recognize the separate ques- 
tion of the convergence of either time or ensemble aver- 
ages. Independent of whether the two averages agree, it 
is possible that the nature of athermal fluctuations in a 
driven system are such that average quantities are not 
well-defined on experimentally accessible time scales or 
number of ensembles. This is highlighted by work in 
granular matter in which extremely long times were re- 
quired for convergence of the average density under tap- 
ping fXlj ] . As we will demonstrate, for our system, both 
the time and ensemble averages are well-defined. 

In this Letter, we report on measurements of the aver- 
age velocity profile using a model two-dimensional foam: 
a bubble raft [l2| • One advantage of this system is the ir- 
relevance of thermal fluctuations due to the macroscopic 
nature of the bubbles. Instead, the fluctuations are the 
result of bubble rearrangements during flow fl3j . Even 
though the average velocity profile for foam has been the 
subject of substantial theoretical 3, 15, H 13] and ex- 
perimental 18|, EI H ED, l23j work, the connection 
between time and ensemble averages for velocity profiles 
remains an open question. 

The bubble raft is produced by flowing regulated ni- 
trogen gas through a needle into a homogeneous solution 
of 80% by volume deionized water, 5% by volume Miracle 
Bubbles (from Imperial Toy Corporation), and 15% by 
volume glycerine. The glycerine provides additional sta- 
bilization of the bubble raft as a function of time. The 
bubbles are confined between two parallel bands sepa- 
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rated by a distance d = 57 mm. The bands are driven 
at a constant velocity v w in opposing directions. This 
applies a steady rate of strain to the system given by 
7 = 2v w /d = 1.4 x 10~ 3 s _1 . The total applied strain is 
7 = jt, where t is the time interval under consideration. 
The direction of flow imposed by the bands is taken to 
be parallel to the x-axis. We use a CCD camera to cap- 
ture the state of the system at a frame rate of 10 Hz. A 
Particle Image Velocimetry program 24j is then used to 
extract velocities of individual bubbles using consecutive 
images. Each consecutive pair of images corresponds to 
a measurement of its instantaneous velocity. For this pa- 
per, we focus exclusively on the behavior of the velocity 
in the x-direction, and the system was divided into forty 
bins in the y-direction so as to measure the velocity as a 
function of y. The bubbles are relatively monodisperse, 
with a bubble diameter of D = 2.66 ± 0.2 mm. Details 
of the apparatus and methods can be found in Ref. [25| . 

The results presented in this letter are based on thirty 
different experimental realizations of a bubble raft. All 
of the parameters for each experimental realization were 
the same except for the configuration of the bubbles. In 
each of the experimental realizations, the system is sub- 
jected to a total strain of five and N = 1000 consecutive 
images of the bubbles are captured. From these thirty 
different experimental realizations, we develop two dif- 
ferent data sets: a time set and an ensemble set. The 
time set is the 30 independent time series of the bubble 
velocities. A member of the time set corresponds to a 
single experimental realization with 1000 elements. Each 
element captures the instantaneous position and veloc- 
ity of the bubbles as measured between successive im- 
ages at a given time (or value of strain). A member of 
the "ensemble set" is constructed from the thirty experi- 
mental realizations by randomly sampling 1000 elements 
from the thirty experimental realizations. We build 30 
such members to constitute the ensemble set. This pro- 
cedure ensures an equal number of elements in each of 
the members of the time and ensemble sets. This allows 
for a direct statistical comparison between the two data 
sets. A time/ensemble average refers to averaging all in- 
stantaneous velocities from a individual member of the 
time/ensemble set of data. 

To characterize the system's behavior, we measured 
the distribution of velocities, the time-correlation of the 
velocities, the average velocity profiles, and characterized 
the convergence rate of the averaging process. All of 
these studies (except the velocity correlations) were made 
for both the time and ensemble data sets. The results for 
the velocity distributions are shown in Fig. [T] for a single 
member of both the time and ensemble set. Sampling 
the velocities from either time sets (26| or ensemble sets 
result in a Lorentzian distribution. However, as we will 
show, a key distinction lies in the distribution of the mean 
values for different members of the respective sets. 

The average velocity profiles are approximately linear 
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FIG. 1: (color online) The probability distribution of veloc- 
ities is indicated for ensemble averaging (open circles) and 
time averaging (closed squares). The solid lines are Lorenz- 
tian fits to the data, red curves for time averaged data and 
blue for ensemble averaged data, (a) Results for a position at 
the center of the trough. For the ensemble (time) average the 
center of the Lorentzian fits are at 0.01 ± 0.01 (-0.01 ± 0.01) 
and the widths are 0.68 ± 0.02 (0.72 ± 0.02), respectively, (b) 
For a position 0.25d from the center. For the ensemble (time) 
average the center of the Lorentzian fits are at 0.41 ± 0.01 
(0.36 ± 0.01) and the widths are 0.47 ± 0.03 (0.52 ± 0.02), 
respectively. 



[25| . To characterize the relative deviations of the various 
velocity profiles in different sets, Fig. [2] plots the differ- 
ence between the individual average velocity profiles and 
the mean profile (v m ) for ten members of the time set 
and ten members of the ensemble set. In the case of time 
averaging, the resulting average velocity profiles qualita- 
tively differ from each other and exhibit relatively (com- 
pared to the ensemble profiles) large deviations from the 
mean. In contrast, the ensemble averages are essentially 
indistinguishable from each other (this is quantified by 
the standard deviation between profiles, see Fig. [2]). 

There are two key issues regarding the average pro- 
files: (1) are they computed over independent configu- 
rations? and (2) have the averages converged? Fig. [3] 
and [4] address these two questions. In Fig. [3l the au- 
tocorrelation of instantaneous velocities for a member 
of a time set and an exponential fit to the correlation 
function (oc exp(— 7/O.OI)) are shown. The decay con- 
stant corresponding to a strain of 0.01 is consistent with 
the previously reported yield strain in polydisperse bub- 
ble rafts (~ 0.05) [27| • As the yield strain corresponds 
to the point at which elastic response gives way to flow 
through significant bubble rearrangements, it is reason- 
able that this determines the interval of strain necessary 
to produce independent configurations. Therefore, time 
averages using a total strain of 5 consist of between 100 
to 500 independent configurations. While building an 
ensemble member, we randomly choose states of veloc- 
ity configurations across all thirty members of the time 
set. This results in the average strain between randomly 



3 




-10 -5 5 -5 5 10 



y/D y/D 

FIG. 2: (a) Deviation of time average velocity profiles (nor- 
malized by the wall speed v m ) from the mean profile (v m ) 
generated by ten members of the time set. (b) Deviation 
of ten different ensemble averages of the velocities from the 
mean profile. Plotting the deviation highlights the differences 
between the different time average profiles. The average stan- 
dard deviation for time average profiles in (a) is 0.081, in 
contrast, the ensemble averages shown in (b) are relatively 
indistinguishable with a standard deviation of 0.017. The y- 
position is scaled by the bubble diameter D. 
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FIG. 3: This is the decay in the correlation of the mean ve- 
locity 0.25d away from the trough center. The total strain 
in each experiment is 5, and we only show the autocorrela- 
tion to a unit strain to make it visually easier to read. The 
remainder of the curve is qualitatively unchanged at higher 
strains. The inset shows a closeup at short strain along with 
an exponential fit (correlation oc exp(— 7/O.OI). 

chosen states being 0.15, with a correlation less than the 
baseline noise in Fig. [3] The selection of these states 
across 30 independent time members, further decreases 
any correlations between them. 

We find the rate of convergence to the mean value 
of the velocity profile is the same for members of the 
time and ensemble sets. However, the distribution of 
final mean values attained are different, with the time 
set showing a significantly larger spread in the converged 
mean values than those from the ensemble set. In Fig. 2J 
the convergence of the mean velocity as a function of the 
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FIG. 4: (color online) The mean values of the velocity 0.25d 
from the center as a function of the number of bubbles used 
to compute the average for (a) time averaging and (b) en- 
semble averaging. For illustration purposes, 10 of the 30 runs 
are shown. The deviation in the final values (indicated by 
the arrows and 5) is larger for time averages than ensemble 
averages. When we use the complete set of 30 runs, we ob- 
tain 8 = 0.258mm/s for time averages and 5 = 0.052mm/s for 
ensemble averages. 

number of bubbles used to compute the average is plot- 
ted for ten members of the time and ensemble averages. 
Plotting the mean velocity versus the number of bubbles 
allows for direct comparison of the time and ensemble 
averages. However, it should be noted, that for the time 
averages, the number of bubbles corresponds directly to 
increasing time, or strain, which in this case is a strain 
< 5. In all cases (ensemble and time averages), the mean 
values converge to a well-defined value after averaging 
over approximately n ~ 2000 bubbles. (This corresponds 
to a strain of 1 for the time averages, or approximately 
100 independent configurations, and is consistent with 
previous measurements of the rate of convergence [l9|). 
The spread in the final mean values is very different for 
the two sets. For the time averages, as expected from 
Fig. [2J the deviation in final mean values is larger than 
the fluctuation in the mean for large enough values of n. 
For the ensemble averages, the deviation in final mean 
values is similar to the fluctuation in the mean. To quan- 
tify this deviation, for a given set of average velocities, we 
define the deviation 8 = v max - v min , where v max (v min ) 
is the largest (smallest) final average velocity in the set. 

Next, in Fig. [5] we describe the spread in the final mean 
values of the velocity profile parameterized by S as a func- 
tion of an increasing number of independent members of 
the time set, M, used in constructing the ensemble set. 
This is done to test whether or not thirty independent 
members of the time set is sufficient to build an ensem- 
ble set with physical relevance, i.e. uncorrelated and in- 
dependent. One method of testing this is to consider 
ensembles constructed from constructed M < 30 mem- 
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FIG. 5: (color online) Plot of the average deviation < S(M) > 
between different ensemble averages as a function of M, the 
number of runs used to generate the ensemble set. As the 
number of independent runs used to generate the ensemble in- 
creases, the deviation between computed averages decreases. 



bers of the time set, and measure S(M). To increase the 
statistics, a number of ensemble sets are generated for 
each M < 30, and one considers < 5{M) > over these 
sets. Figure[5]is a plot < S(M) > for each M, and it con- 
firms our expectation that increasing M should increase 
the independence between the randomly selected veloci- 
ties, resulting in a better approximation of an ensemble 
average. For M > 20, the value of < S(M) > is bounded 
within range of the error bar of the value at M = 20. This 
has two important consequences. First, the independence 
of < 6(M) > on M for M > 20 is strong evidence that 
thirty independent experimental realizations provides a 
sufficient approximation of the ensemble average for this 
system. Second, the indication that < S(M) > has con- 
verged to a nonzero value may be attributed to the finite 
size of the system. An interesting future study will be 
the dependence of < S(M) > as a function of the system 
size. 

We have presented three main results. First, ather- 
mal fluctuations generated during the flow of a bubble 
raft are sufficient to produce stationary, time-averaged 
velocity profiles after a total strain of approximately one 
(Fig. 2]). Second, these fluctuations do not provide suf- 
ficient exploration of phase space to result in a unique 
time-average. In other words, different realizations of a 
bubble raft produce different time-averaged velocity pro- 
files that exhibit a measurable and significant deviation 
from each other (see Fig. 2] and [5]) . Third, by sampling 
velocities from different realizations of the experiment, 
the resulting ensemble averages demonstrate a conver- 
gence towards similar velocity profiles (Fig. [4] and [2]). 
Putting these results together, one concludes that the 
limit of time averaging is different from that of ensemble 
averaging for this system in at least one crucial aspect 
i.e., there is no convergence to a unique time-average for 
all experimental realizations of the bubble raft, but the 
ensemble averages converge to final velocity profiles that 
are consistent with a unique profile. 
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